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The detection of the theoretically expected dark matter is central to particle physics 
and cosmology. Current fashionable supersymmetric models provide a natural 
dark matter candidate which is the lightest supersymmetric particle (LSP). The 
allowed parameter space of such models combined with fairly well understood 
physics (quark substructure of the nucleon and nuclear structure) permit the eval- 
uation of the event rate for LSP-nucleus elastic scattering. The thus obtained 
event rates, which sensitively depend on the allowed parameter space parameters, 
are usually very low or even undetectable. So, for background reduction, one would 
like to exploit two nice features of the reaction, the directional rate, which depends 
on the sun's direction of motion and the modulation effect, i.e. the dependence 
of the event rate on the earth's annual motion. In the present paper we study 
these phenomena in a specific class of non isothermal models, which take into ac- 
count the late in-fall of dark matter into our galaxy, producing flows of caustic 
rings. We find that the modulation effect arising from such models is smaller than 
that found previously with isothermal symmetric velocity distributions and much 
smaller compared to that obtained using a realistic asymmetric distribution with 
enhanced dispersion in the galactocentric direction. 



I. Introduction 

In recent years the consideration pf exotic dark matter has become necessary 
in order to close the Universe a> 0. Furthermore in in order to understand 
the large scale structure of the universe it has become necessary to consider 
matter made up of particles which were non-relativistic at the time of freeze 
out. This is the cold dark matter n component (CDM). The COBE data B 
suggest that CDM is at least 60% On the other hand during the last few 
years evidence has appeared from two different teams, the High-z Supernova 
Search TeamB and the Supernova Cosmology Project □ ' El which suggests that 
the Universe may be dominated by the cosmological constant Ar. As a matter 
of fact recent data the situation can be adequately described by a barionic 
component tts = 0.1 along with the. exotic components CIcdm — 0.3 and 
f?A = 0.6. In another analysis TurnerO gives Vt m = tie dm i+ &b = 0.4. Since 
the non exotic component cannot exceed 40% of the CDM H'ta, there is room for 
the exotic W-IMP's (Weakly Interacting Massive Particles). In fact the DAM A 
experiment 113 has claimed the observation of one signal in direct detection of 
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a WIMP, which with better statistics has subsequently been interpreted as a 
modulation signal til. 

The above developments are in line with particle physics considerations. 
Thus, in the currently favored supersymmetric (SUSY) extensions of the stan- 
dard model, the most natural WIMP candidate is the LSP, i.e. the lightest 
supersymmetric particle. In the most favored scenarios the LSP can be simply 
described as a Majorana fermion, a hii£ar_combination of the neutral compo- 
nents of the gauginos and Higgsinoscrlid til. 

Since this particle is expected to be very massive, m x > 30GeU, and 
extremely non relatipstic with average kinetic energy T < lOOKeV, it can be 
directly detected t3~lij mainly via the recoiling of a nucleus (A,Z) in the elastic 
scattering process: 



X + (A,Z) -> x + (Azy (1) 

(x denotes the LSP). In order to compute the event rate needs the following 
ingredients: 

1) An effective Lagrangian at the quark-level ip-ithe context of supersym- 
metry as described in Rcfs.u, Bottino et al. Iii3 and til. 

2) A a quark model for the nucleon, needed in going from the quark to the 
nucleon level, since the-ebtained results are sensitive to the presence of quarks 
other than u and d liilaH 

3) Compute the differentiaLcioss sections using as reliable as possible many 
body nuclear wave functions. EirEZl. 

The obtained rates sensitively depend on the input from the allowed SUSY 
parameter space. Since the expected rates are extremely low or even unde- 
tectable with present techniques, one would like to exploit the characteristic 
signatures provided by the reaction. Such are: a) The modulation effect, i.e 
the dependence of the event rate on the velocity of the Earth and b) The 
directional event rate, which depends on the velocity of the sun around the 
galaxy as welLas the the velocity of the Earth. The latter effect, recognized 
sometime agac3 has recently begun to appear feasible by the planned UKDMC 
experiment Ej. We will study both of these effects in the present work. 

In our previous letter E^land its subsequent expanded version 113 we found 
enhanced modulation, if one uses an appropriate asymmetric velocity distri- 
bution with enhanced dispersion in the galactocentric directionEU. 

The isolated galaxies are, however, surrounded by cold dark matter—which, 
due to gravity, keeps falling continuously on them from all directions E3. As a 
result one has caustic rings with matter density, which depends on space and 
velocity. 

It is the purpose of our present paper to exploit the results of such a sce- 
nario and calculate quantitavely the resulting modulation effect in the usual 
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(differential and total ) rate for LSP-nuclcus clastic scattering. We will also 
study the directional rates in the spirit of the recent work of Copi et alE3. 
For the reader's convenience we will give a very brief description of the basic 
ingredients on how to calculate LSP-nucleus scattering cross section. We will 
not, however, elaborate on how one gets the needed parameters from super- 
symmetry. The calculation of these parameters has become pretty standard. 
One starts with representative input in the restricted SUSY parameter space 
as described in the literature, e.g. Bottino et al. 113, Kane et al. , Castano et 
al. and Arnowitt et al. I1!j. 

After this we will specialize our study jn .the case of the nucleus 127 1, 
which is one of the most popular targets Efr . We will present our results 



a function of the LSP mass, m x , in a way which can be easily understood by 
the experimentalists. 



II. The Basic Ingredients for LSP Nucleus Scattering 

Because of lack of space we are not going to elaborate here further on the 
construction of the effective LaffiraugLan derived from supersymmetry, but refer 
the reader to the literatuyEj'y^EjilatJ. For the reader's convenience we will 
summarize our previous E-3 formulas related to the LSP-Nucleus cross section 
and the event rates. The effective Lagrangian can be obtained in first order 
via Higgs exchange, s-quark exchange and Z-exchange. In a formalism familiar 
from the theory of weak interactions we write 

Leff = -^|{(xi7 A 75Xi)Ja + (xiXi)J} (2) 

where 



Jx = N lx {f v + f v r 3 + f Al 5 + /A7 5 r 3 )iV (3) 

and 

J = N(fs + fln)N (4) 

We have neglected the uninteresting pseudoscalar and tensor currents. 
Note that, due to the Majorana nature of the LSP, Xi7 A Xi = (identically). 
The parameters f v , f v , f\, f\, fg, fg depend on the SUSY model employed. 
In SUSY models derived from minimal SUGRA the allowed parameter space 
is characterized at the GUT scale by five parameters, two universal mass pa- 
rameters, one for the scalars, mo, and one for the fermions, nil /?, as well as 

the parameters tan/3, one of Aq (or m\ ole ) and the sign of /xE3. Deviations 
from universality at the GUT scale have also been considered and found useful 
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0. Wc will not elaborate further on this point since the above paramcUp. 
involving universal masses have already been computed in some models Mr 6 ! 
and effects resulting from deviations from universality.-will be found elsewhere 
(see Arnowitt et al in Ref . cJ and Bottino et al in Ref. Ii3) . For some choices in 
the allowed parameter space the obtained couplings can be found in a previous 
paperE3. 

With the above ingredients the differential cross section can be cast in the 
form 



du 



Mu,v) = ^ _ [(Sg + S y ~zrz5nz~) F 2 (u) + t spm F 11 {u)} (5) 



E S = a (J^r {A 2 [(f°s - ti^-^f } - al xa A\^f-f (6) 



y\ Spin J- /ry\ 

2->spin — 0"p ^-0 <,spin \ * } 

Cs Pm = ( ^^y ))2 [(f» (o)) 2 g^ + 2^O (0)O l( 0)fiM + n l( o))»] 

3(1 + 4)2 fX F n {u) fX F n {u) 



(8) 

t v = a^ x o (v (9) 



_ (Vr/MN)) 2 ,a M A-2Z 2 <t; 2 > 1 2t ? +1(2 M ) . 

^ (1 + 4) 2 /{/ A c 2 (2^ r b) 2 (l + V ) 2 (v 2 )- 



(10) 



a* %0 = proton cross-section, i = 5, spin, V given by: 

a px° = a ° (fs) 2 ( sca l ar ) > (the isovector scalar is negligible, i.e. = <r„) 



a 



p,x" 

spin 



,o - to 3 (/° + /i) 2 (spin) , a* = ao (/° + fy) 2 (vector) 



where m p is the proton mass, rj = m x /m.NA, [i r is the LSP-nucleus reduced 
mass , [i r (N) is the LSP-nucleon reduced mass and 

^0 = ^-(G F m N ) 2 ~ 0.77 x l(T 38 cm 2 (11) 

Z7T 



with 



u = q 2 b 2 /2 or Q = Q u, Q = _L_ (12) 

where b is (the harmonic oscillator) size parameter, q is the momentum transfer 

and Q is the energy transfer to the nucleus 

In the above expressions F (u) is the nuclear form factor and 

are the spin form factors 111 (p, p are isospin indices) Both form factors are 
normalized to one at u = 0. flo (^1) are the static isoscalar (isovector) spin 
matrix elements. 

The non-directional event rate is given by: 

dN p(0) m . .. 

R = Rnon-dir = —jr = 1 a(u,V)\v\ (14) 

at m x ArriN 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity and m is the 
detector mass 

The differential non-directional rate can be written as 

dR = dR non ^ dlr = da(u,v)\v\ (15) 

m x ArriN 

where do~(u, v) was given above. 

The directional differential rateli!3 in the direction e is given by : 

dRdir = . m v.eHOu.e) ^- da(u,v) (16) 

where H the Heaviside step function. The factor of 1/27T is introduced, since 
the differential cross section of the last equation is the same with that entering 
the non-directional rate, i.e. after an integration over the azimuthal angle 
around the nuclear momentum has been performed. In other words, crudely 
speaking, l/(2ir) is the suppression factor we expect in the directional rate 
compared to the usual one. The precise suppression factor depends, of course, 
on the direction of observation. 
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III. Convolution of the Event Rate 

We have seen that the event rate for LSP-nucleus scattering depends on the 
relative LSP-target velocity. In this section we will examine the consequences 
of the earth's revolution around the sun (the effect of its rotation around its 
axis is expected to be negligible) i.e. the modulation effect. In the past this 
has been accomplished by .assuming a consistent LSP velocity dispersion, such 
as a Maxwell distribution □ or asymmetric jaelocity distribution with enhanced 
dispersion in the galactocentric direction Eiliilltll. More recently other very 
interesting non-isothermal approaches have been proposed, which consider the 
in- fall of dark matter into the galaxy producing flows of caustic rings. In 
particular the predictions of a self-similar model have been put forward as a 
possible scenario for dark matter density-velocity distribution, see e.g. Sikivic 
et all£3. The implications of such distributions on the direct detection rates 
and, in particular, on the directional ones and the modulation effect are the 
subject of this work. Before proceeding furtyher we should mention that after 
our manuscript had been prepared another approach has been suggested, which 
links the distribution to the density profile via Eddigton's formulaEl. 

Following Sikivie we will consider 2 x N caustic rings, (i,n) , i=(+.-) and 
n=l,2,...N (N=20 in the model of Sikivie et al), each of which contributes to 
the local density a fraction p n of the of the summed density p of each of the 
£ = +,—. It contains WIMP like particles with velocity y n = {y nx ,V ny ,y nz ) 
in units of essentially the sun's velocity (vq — 220 Km/ s), with respect to the 
galactic center. The z-axis is chosen in the direction of the disc's rotation, i.e. 
in the direction of the motion of the the sun, the y-axis is perpendicular to the 
plane of the galaxy and the x-axis is in the radial direction. We caution the 
reader that these axes are traditionally indicated by astronomers as e^, e r , e^X 
respectively. The needed quantities are taken from the work of Sikivie et alE3 
(see in TABLE I), via the definitions 

V'n = v n/v ,y nz = V ntj ,/v a = y n z,y'nx = v nr/v = Vnx^Vny = V nz /v = 

y ny ,Pn = d n /p, p = J2n=i d n and Vn = [{Vnz - I) 2 + vl y + vl x } 1/2 (for each 
flow +.-). This leads to a velocity distribution of the form: 

N 

f(v>) = £ S(v - v y' n ) (17) 

n=l 

Since the axis of the ecliptic El lies very close to the y, z plane the velocity of 
the earth around the sun is given by 

Ve = Vo + t>i = wo + ui ( sina x — cosa cosj y + cosa siwy z ) (18) 

where a is the phase of the earth's orbital motion, a = 2n(t — t±)/TE, where 
t\ is around second of June and Tg = lyear. 
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One can now express the above distribution in the laboratory frame 0by 
writing v = v + ve 

IV. Expressions for the Differential Event Rate. 

The mean value of the non-directional event rate of Eq. (|l5|), is given by 
/dR\ p(0) m f da(u,v) l3 , 

(*T> = ifl^ J KwM-jb 1 ** (19) 

The above expression can be more conveniently written as 



where 



,dS x f \v\ tl s dcr{u,v) 3 



V ^ (21) 



c?w 7 \/(v 2 ) ' 

There are now experiments under way aiming at measuring directional 
rates using TPC counters which permit the observation of the recoiling 
nucleus is observed in a certain direction. From a theoretical-point of view the 
directional rates have been previously discussed by SpergelEj and Copi et al 
EJ. The mean value of the directional differential event rate of Eq. ([l6|), is 
defined by 

/dR\ p(0) m 1 f ,, . „„, .d<j{u,v) -> . . 

( — ) =^- L - A / f(v,v E )v.eH(v.e) y , ' d 3 v (22) 

\ du I dir m x ArriN 27r J du 

where e is the unit vector in the direction of observation. It can also be more 
conveniently expressed as 



\ du I dir m x ArriN du dv 



where 



,cK. 1 f v.eH(v.e) „, .da(u.v) ,n 

< = J ^7PT /(V ' VE) ^dt ld v (24) 

It is clear that the rate will depend on the direction of observation, showing 
a strong correlation with the direction of the sun's motion. In a favorable 
situation the rate will merely be suppressed by a factor of 2tt relative to the 
non-directional rate. We will specialize our results in the case of caustic rings. 
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A. The Non- directional Differential Event Rate in the Case of caustic Rings 
Eq. El] takes the form 



, dS 2p 



du' p(0) 



a'psfoW + ^SyFi(w) + Y> apin F spin (u)] (25) 



We remind the reader that p was obtain for each, type of flow (+ or -), which 
explains the factor of two. In the Sikivie modelE3 we have (2p/p(0) ~ 1.0, i.e. 
the whole dark matter density lies in the form of caustic rings. In hybrid mod- 
els, which contain in addition an isothermal component, it is only a fraction, 
since the sum of all densities should be p(0). 

The quantities E$,i = S,V,spin are given by Eqs. (^|)- (0). The quan- 
tities F@ , F\ , F S pi n are obtained from the corresponding form factors via the 
equations 

F k (u)=F 2 (u)* k (u)Q±^ , fc = 0,l (26) 



F spin {u) =F n (u)y (u) (27) 
The functions ^k{u) depend on the model. Introducing the parameter 

S = — = 0.27, (28) 

we find in the case of Sikivie model 

AT 



n=l 

with 



^/2p r bvo 

and 

xi & ■ \1 i & \2 , i & 

Vn = [{Vnz ~ 1 - ■x sin^i coset) + (y ny + - cosj cosa) + (y n 



— stna 



(29) 



(30) 



.211/2 



2 / / ' \a-ay T 2 * ^ IJC "- 1 / ~ \ynx ^ 

(31) 

The above expressions (p5|)-(p7|) depend not only on u, but on the phase of 
the Earth as well as in the angle 7 ~ ir/2. If, however, in the 8-function we 
approxmate y n by y n given by 



yn = [(y nz -i)) 2 + y 2 ny + yL} 1/2 (32) 

i.e. neglect threshold effects originating from the motion of the Earth, the 
obtained expressions can be simplified. In fact, if we then expand ^(w) in the 
small parameter 6 and keep terms up to linear in it, the dependence on the 
phase of the Earth can be seperated from the dependence on u. Thus, to an 
approximation which is valid to no worse than 10% for all u, the non-directional 
differential rate takes the form 

A = R^f-t T(u)[l- cos a H(u)] (33) 
du p[0) 

where now T(u) and H(u) are now function of the variable u. 

The last expression must be compared to that of isothermal models: 

,dR, - p (0) , w , ., , 

— )=RiL±±t T (u 1 + cosa H (u) 34 
du p{uj 

with p (0) being the density associated with the isothermal mode. In hybrid 
models it is less than p(0).,-Jp the above expressions R is the rate obtained in 
the conventional approachl!3 by neglecting the folding with the LSP velocity 
and the momentum transfer dependence of the differential cross section, i.e. 
by 

*=^^V^Ps + W^] 05) 

where Sj,i = S, V, spin have been defined above, see Eqs (||) - @. 

The factor T(u) takes care of the u-dependence of the unmodulated dif- 
ferential rate. It is defined so that 

duT(u) = 1. (36) 

i.e. it is the relative differential rate. u m i„ is determined by the energy cutoff 
due to the performance of the detector, i.e 

Qmin ('Q7 1 \ 

while u max is determined the via the relations: 

2 2 

u m a X = mm(^f,max(^) , n = 1, 2, N) (38) 
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On the other hand H (u) gives the energy transfer dependent modulation am- 
plitude (relative to the unmodulated amplitude). The quantity t takes care 
of the modification of the total rate due to the nuclear form factor and the 
folding with the LSP velocity distribution. Since the functions F${u),F\ and 
F S pi n have a different dependence on u, the functions T(u)and H (u) and t, in 
principle, depend somewhat on the SUSY parameters. If, however, we ignore 
the small vector contribution and assume (i) the scalar and axial (spiftX de- 
pendence on u is the same, as seems to be the case for light systems L3t3, or 
(ii) only one mechanism (S, V, spin) dominates, the parameter R contains the 
dependence on all SUSY parameters. The parameters t and T(u) depend on 
the LSP mass and the nuclear parameters, while the H(u) depends only on 



the parameter a. Eq. 21 takes the form 



B. Expressions for the Directional Differential Event Rate 

The model of Sikivie et alEU, which is not spherically symmetric, offers itself 
as a perfect example fot-the study of the directional rates. Working as in the 
previous section we get 113 

,dT, 2p a 2 - (v 2 ) - - 

(-j-)dir = tttt tH^s-FoW H 5-Sv-Fi(«) + T, spin F spin {u)\ (39) 

au p[0) 2tt c z 

where the Ej,z = S,V,spin are given by Eqs. (||)- (||). The quantities 
F , Fi , F spin are obtained from the equations 

F k (u) = F 2 (u)* k (u) { j±^-,k = 0, 1 (40) 



F spin {u) =F u (u)*o(u) (41) 
In the Sikivie model we find 

[2 N v 2 S 

= Y g P»Vn k 1)Q (^ ~ u )\(Vnz ~ 1 ~ 2 Sinl C0Sa ^ ez - e 

S 6 
+ (Vny + 2 cos T cosa)e v .e + (y nx - - sina)e x .e\ (42) 

In the model considered here the z-component of the LSP velocity, with respect 
to the galactic center, for some rings is smaller than the sun's velocity, while 
for some others it is larger. The components in the y and the x directions are 
opposite for the + and - flows. So we will distinguish the following cases: a) e 
has a component in the sun's direction of motion, i.e. < 9 < ir/2, labeled by 
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u (up), b) Detection in the direction specified by tt/2 < 9 < n, labeled by d 
(down). The differential directional rate takes a different form depending on 
which quadrant the observation is made. Thus, keeping terms up to linear in 
S, we find: 

1. In the first quadrant (azimuthal angle < <p < 7r /2). 

( ^7 } = R W)^ Tiu)m{u) ~ c ° sa Hi{u))lezM 

+ (R l y + cosaH l 2 (u) H ^-(\cosa\ + cosa))\e y .e\ 

+ (R x - sinaHi(u) H y^(|sma| - sina))\e x .e\] (43) 

2. In the second quadrant (azimuthal angle tt/2 < <f) < tt) 

fjf>i -On t 

W = R W)^ T{u)[{K{u) ~ c ° sa H{[u)) ^ M 

(u) 

+ (R l y + cosaH^u) H ^-(\cosa\ - cosa))\e y .e\ 

+ (Rl + sinaHi(u) H y^(|sma| + sina))\e x .e\] (44) 

3. In the third quadrant (azimuthal angle ir < <f) < 3n/2). 

^ = R W)^ T{u)[mu) cosa Hl{u))lezM 

+ (R l y - cosaH^u) H ^—-(\cosa\ - cosa))|e y .e| 

+ (Rl + sinaHi(u) H ^-(\sina\ + sina))\e x .e\] (45) 

4. In the fourth quadrant (azimuthal angle 3n/2 < cf) < 2n) 

<^> - ^j)^™^^) - cos « Hl(u))\e z .e\ 



(u) 

+ (R y — cosaH^u) H ^ — (|cosa| — cosa))\e y .e\ 



+ (R x - sinaHl(u) + Hs ^\ \sina\ - sina))\e x .e\\ (46) 

where i — u,d 

By the reasoning given above, if one mechanism is dominant, the parame- 
ters R x , R y , R z , Hi, H 2 , H 3 , H c , H s for both directions u and d depend only on 
[i r and a. They are all independent of the other SUSY parameters. 
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V. The Total Event Rates 



We will distinguish two possibilities, namely the directional and the non direc- 
tional case. Integrating Eq. ( |33| ) in the case of caustic rings we obtain for the 
total non directional rate 

R = Rt^L[^-h{a,Q min )cosa)} (47) 

p(o) 

to be compared with the corresponding one for isothermal models: 

R = Rt^[l + h(a,Q mm )cosa)} (48) 

In the above expressions Q m in is the energy transfer cutoff imposed by the 
detector. The modulation is described by the parameter h only. 

The effect of folding with LSP velocity on the total rate is taken into 
account via the quantity t. The SUSY parameters have been absorbed in 
R. From our discussion in the case of differential rate it is clear that strictly 
speaking the quantities t and h also depend on the SUSY parameters. They 
do not depend on them, however, if one considers the scalar, spin etc. modes 
separately. 

Let us now examine the directional rate. Integrating Eqs. (|43| ) - (|6|) we 
obtain: 

1. In the first quadrant (azimuthal angle < <f> < tt/2). 
Rdir = ~ C0Sa h l)\ e - e \ 

h* 

+ (r y + cosah l 2 H — ^(|cosa| + cosa))\e y .e| 

2. In the second quadrant (azimuthal angle 7r/2 < (f> < tt) 

R^ r = R^r^[{r i z -co S ahi)\e z .e\ 

+ (r l y + cosah l 2 (u) + -^-(\cosa\ — cosa))\e y .e\ 

+ (r* + sinah\ + -^-{\sina\ + sma))|e x .e|] (50) 

3. In the third quadrant (azimuthal angle 7r < <fi < 3tt/2). 
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R, 



'dir 



= R 



2p t 
p(0) 2^ 



[{r\ — cos a ft, l 1 )|e z .e| 



+ (r l y — cosah 2 (u) H — ^ — (|cosa| — cosa))\e y .e\ 
+ (r x + sinaH^ + -^-(\sina\ + sina))\e x .e\] 




(51) 



4. In the fourth quadrant (azimuthal angle 3ir/2 < <fi < 2ir) 



2p t 
p(0) 2^ 



{(r l z — cos a /i^)|e z .e| 



+ (rl — cosah l 2 + -^(\cosa\ 
y 2 

h* 

+ (r l x — sinah\ + -^-(\sina\ 



cosa))\e y .e| 



sina))\e x .e 



I] 



(52) 



VI. Discussion of the Results 

We have calculated the differential as well as the total event rates (directional 
and non directional) for elastic LSP-nucleus scattering using realistic nuclear 
form factors. We focused our attention on those aspects of the problem, which 
do not depend on the parameters of supersymmetry other than the LSP mass. 
The parameter R, normally calculated in SUSY theories, was uat calculated 
in this work. The interested readep-4s petered to the literature La ' E3 and, in 
our notation, to our previous worktHl 'EJ E3. 

We specialized our results for the target 127 1. We considered the effects 
of the detector energy cutoffs, by studying two typical cases Q m in — 10, 20 
KeV. Only the coherent mode due to the scalar interaction was considered. 
The spin contribution will appear elsewhere. 

Special attention was paid to the modulation effect due to the annual 
motion of the earth. We assumed that the LSP density in our vicinity and the 
velocity spectrum is that of caustic rings resulting from the self- similar model 
of Sikivie et alE3 . 

We will primarily concentrate on the total rates, which are described 
in terms of the quantities i, r^,,r*,r* for the unmodulated amplitude and 
h, h\, h l 2 , h\, h l c , h l s for the modulated one. In TABLE II we show how these 
quantities vary with the detector energy cutoff and the LSP mass. Of the 
above list only the quantities t and h enter the non directional rate. We notice 
that the usual modulation amplitude h is small. The main reason is that there 
are cancelations among the various rings, since some rings are characterized 
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by Unz > 1) while for some others y nz < 1 (see TABLE I). Such cancelations 
are less pronounced in the isothermal models. As expected, the parameter 
t, which contains the effect of the nuclear form factor and the LSP velocity 
dependence, decreases as the reduced mass increases. 

We observe that the quantities h, ,hj, i = u,d and j = x,y,z,c,s, being 
the ratio of two amplitudes, are here essentially independent of the energy cut- 
off Qmin- On the other hand the quantity t decreases after Q m in is introduced, 
since an important part of the phase space is excluded. 

We notice that, unlike the isothermal models, the maximum of modulation 
oqeurs around the 2nd of December, something already noticed by Sikivie et 
aM 

Let us now examine the differential rates. We will begin with non-directional 
one, which is described in terms of the functions T(u) and H(u). These are 
shown for various LSP masses and Qmin in Fig- la for T(u) and Fig. lb for 
H(u). We remind the reader that the dimensionless quantity u is related to 
the energy transfer Q via Eq. ( 113 ) with Qq — 60KeV for 127 1. Note that, due 
to our normalization of T, the area under the corresponding curve is unity. 
This normalization was adopted to bring the various graphs on scale, since 
the absolute values may vary substantially as a function of the reduced mass. 
We observe that the function T(u) differs somewhat from the predictions of 
the isothermal models models. Here the function begins with a maximum at 
u = 0.0, while in the isothermal models the maximum occurs at u = 0.1. Fur- 
thermore this function shows less of a reduction as the reduced mass increases 
(see Fig la). The functions H(u) exhibit step behavior in some regions. This 
is not unexpected, since the LSP velocity spectrum was assumed to be discreet 
in the Sikivie model. 

The directional differential rates, which is now beginning to look like a 
feasible experimental possibility E3, is described by the functions 

R x ,R y , R z , Hi,H 2 , H 3 , H c , H s 
(for both directions u and d) Due to lack of space we are not going to show the 
here. We only mention that they were defined as ratios of functions, with T{u) 
in the denominator. With this definition they are essentially independent of 
u, with values approximately equal to their corresponding total values shown 
in TABLE II. 

VII. Conclusions 

In the present paper we have calculated the parameters, which describe the 
event rates for direct detection of supersymmetric dark matter. We studied, 
in particular, the directional variation of the rates and the modulation effect. 

The needed local density-iand velocity spectrum of the LSP were taken 
from the work of Sikivie et alE3, viewed as a late in-fall of dark matter into our 



14 



galaxy. They were derived in the context of a self-similar model, which yields 
40 caustic rings. 

We presented our results in a suitable fashion, i.e by separating the rates 
into two factors. One factor R, which carries the dependence on the SUSY 
parameters , not the subject of the present work, and another, which is essen- 
tially independent of all SUSY parameters except the LSP mass. The latter 
depends mainly on the properties of the LSP velocity distribution, the nu- 
clear parameters and the kinematics. The nuclear form factor was taken into 
account and the effects of the detector energy cut off were also considered. 
Strictly speaking the obtained results describe the coherent process in the case 
of 127 1, but we do not expect large changes, if the axial current is considered. 

Our results, in particular the parameters t, see TABLE II, indicate that 
for large reduced mass, the kinematical advantage of /i r (see Eqs. (||)- ([h]) is 
partly lost when the nuclear form factor and the convolution with the velocity 
distribution are taken into account. Also, if one attempts to extract the LSP- 
nucleon cross section from the data, in order to compare with the predictions 
of SUSY models, one must take t into account, since, for large reduced mass, 
t is different from unity. 

In the case of the non-directional total event rates we find that the max- 
imum no longer occurs around June 2nd, but about six months later. The 
difference between the maximum and the minimum is about 4%, a bit smaller 
than that predicted by the symmetric isothermal models. It is, however, sub- 
stantially-Less, than h = 0.46 predicted by asymmetric isothermal velocity dis- 
tribution H-El. 

In the case of the directional rate we found that the rates depend on the 
direction of observation. The biggest rates are obtained, if the observation 
is made close to the direction of the sun's motion. The directional rates are 
suppressed compared to the usual non-directional rates by the factor f rec i = 
k/(2tt). We find that k = r" ~ 0.7, if the observation is made in the sun's 
direction of motion, while k ~ 0.3 in the opposite direction. The modulation 
is a bit larger than in the non-directional case, but the largest value, 8%, 
is not obtained along the sun's direction of motion, but in the x-direction 
(galactocentric direction) . 
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TABLE I.: The velocity parameters y' n = v n /vo, y nz = y' nz = v n4 ,/v , y ny = y' ny = 

Vnz/vo, y nx = y nx = Vnr/vo and y n = [(y nz - l) 2 + y 2 y + ^J 1 / 2 . Also shown are the 

quantities: a„, the caustic rind radii, and p n = ^n/[X^n— 1 ^ n l' ( ^ or t ' lc otner definitions 
see text ). 
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TABLE II. : The quantities t and h entering the total non directional rate in the case of the 
target 53 1 127 for various LSP masses and Q m i„ in KcV. Also shown are the quantities r*., h z - 
i = u,d and j = x, y, z, c, s, entering the directional rate for no energy cutoff. For definitions 
sec text. 
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FIG. 1.: The quantities T(u) and H(u) entering the differential amplitude. Thick solid line 
corresponds to m c hi = 30 GeV the intermediate thickness line to m x = 80 GeV, the fine 
line to m x = 100 GeV. The rest correspond to larger LSP masses and fall on top of each 
other. f 



21 



